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SUMMARY 


PROBLEM 

Investigate  alternate  means  to  the  split-step  algorithm  for  determining  coverage  in  the  troposphere 
at  gigahertz  (GHz)  frequencies. 

RESULTS 

The  theory  of  the  implicit-finite-difference  method  is  summarized  and  the  boundary  conditions 
described.  Examples  indicate  that,  for  sufficiently  high  signal  levels,  the  IFD  method  is  adequate  at 
9.6  GHz.  However,  the  method  did  not  give  correct  results  beyond  35  km  for  the  standard-atmos¬ 
phere  case  in  which  signal  levels  are  very  low. 
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I.  INTRODUCTION 


This  report  describes  an  effort  to  investigate  alternate  means  to  the  split-step  algorithm  for  determin¬ 
ing  radar  coverage  in  the  troposphere  at  gigahertz  (GHz)  frequencies.  It  concerns  mostly  the  results  of 
testing  the  implicit-finite-difference  (IFD)  algorithm  described  by  Lee,  Botseas,  and  Papadakis  [1981]. 
Both  methods  are  solutions  to  the  parabolic  wave  equation  (PE),  for  which  reflections  are  neglected. 

Lee,  et  al.  [1981],  give  an  analysis  of  truncation  error  for  the  implicit-finite-difference  method  and 
conclude  that  it  is  consistent,  is  stable,  and  converges.  S.  T.  McDaniel  [1975]  concluded  that  the  finite- 
difference  method  is  useful  only  at  low  frequencies.  However,  results  were  obtained  by  this  method  at  9.6 
GHz,  except  for  one  case  in  which  signal  levels  are  very  low. 

Section  II  of  this  report  gives  a  summary  of  the  theory  of  the  implicit-finite-difference  method.  Sec¬ 
tion  III  summarizes  the  starting  conditions  used  at  the  beginning  of  the  integration  of  the  parabolic  equa¬ 
tion.  In  section  IV,  three  examples  are  given  of  use  of  the  IFD  method. 

The  characteristics  of  the  tropospheric  models  are  given  as  modified  refractivity,  M,  versus  height 
where  M  =  (n  -  1)  x  106  +  0.157  z,  n  is  the  index  of  refraction,  and  z  is  the  height  in  meters.  Results  are 
given  in  terms  of  path  loss  versus  distance. 


II.  THEORY  OF  THE  IMPLICIT-FINITE-DIFFERENCE  METHOD 

The  parabolic  equation  i*  of  the  form 


du/dr  =  a(k0,  r,  z)  u  +  b{k0)  d2u/dz2 
=  (a  +  bD 2)  m  =  Lu 


where 


a  =  (ik0/2)\n2{r,z)  -  1] 
b  =  i/2k0 
D  =  d/dz 


and  where  k0  is  the  reference  wave  number,  n  is  the  index  of  refraction,  r  is  range,  and  z  is  altitude.  This 
can  be  expressed  in  the  symmetrical  form 


e~  '/2kL  Uml  =  e  1 ,2kL  Um 


(2) 


where  m  is  an  index  with  regard  to  mesh  points  in  the  vertical  direction  and  n  is  an  index  with  regard  to 
mesh  points  in  the  horizontal  direction.  Claerbout  [1970]  indicates  that  use  of  a  symmetrical  form  of 
solution  results  in  considerably  more  accuracy. 

To  first  order  in  the  exponentials,  equation  2  becomes 


[  1  -k(aml  +  bm'D2)/2]Um'  =  [  1  +  k{a%  +  bnmD2) /2]unm. 


(3) 
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This  approximation  inherently  limits  the  mesh  size.  Nevertheless,  in  two  ducting  cases  at  9.6  GHz,  the 
method  gave  accurate  results  with  reasonably  large  mesh  size. 

Following  Lee,  et  al.  [1981],  equation  3  may  be  written  in  matrix  form  as 
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where 


Xm  =  l-ka^/2^[k/(6z)2]bnm^ 
Ym=  \+kanm/2-[kH6z)2]bnm 
Pnm  =  [*/(&)2]&8,. 


The  nonzero  elements  of  the  last  two  vectors  represent  quantities  to  be  determined  by  boundary 
conditions  at  the  top  and  bottom  of  the  troposphere  model.  The  first  mesh  point  of  the  model  was  taken 
to  be  one-half  mesh  unit  above  the  ocean  surface.  The  zero^  mesh  point  was  therefore  one-half  mesh 
unit  below  the  surface.  The  ocean  was  assumed  to  be  a  perfect  conductor.  Therefore,  the  bottom  bound¬ 
ary  condition  was  u"  =  u". 

The  boundary  condition  at  the  top  should  be  that  only  outgoing  waves  are  present.  This  is  modeled  by 
placing  an  absorber  in  the  upper  part  of  the  profile.  Following  Brock  [1978],  the  upper  one-fourth  of  the 
profile  was  modified  by 


n2  *-  n2  -  /0.01  exp  {- [(z -ztop)/(z(o;,/4)]2} 


where  i  is  J-  1  and  the  difference  in  sign  is  the  result  of  the  electromagnetic  convention  of  waves  traveling 
as  exp(/'a><)  compared  to  the  acoustical  convention  of  waves  traveling  as  exp(-/a>/) 

The  value  of  was  determined  from  the  condition 
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a  is  the  radius  of  the  Earth,  and  <5z  is  the  vertical  mesh  size. 

Another  upper-boundary  condition  was  tried  for  which  the  value  was  taken  to  be  its  value  as 
determined  by  ray  tracing.  In  such  cases,  no  absorber  was  used.  No  such  run  was  successful,  however. 

III.  STARTING  SOLUTION 

The  parabolic  equation  method  requires  a  set  of  starting  values  as  a  function  of  height.  Two  methods 
were  employed.  One  was  to  use  a  “free  space”  condition  except  that  the  Earth  as  a  flat  perfect  conductor 
was  included.  It  was  taken  to  be  valid  only  somewhat  close  to  the  transmitter  and  in  practice  was  used  at 
2.5  km  to  10  km  from  it.  It  is  given  by 

u  =  [exp(ai)  -  exp(a2)]/ J7o 

where 

c*i  =  -ik{[rl  +  (z-zr)2]1/2-r0} 

a  2  =  -ik{[rl  +  (z  +  zr)2]1/2  -r0) 

zf  is  the  height  of  the  transmitter,  and  ra  is  the  range. 

The  other  method  was  a  ray-trace  formulation  developed  by  R.  A.  Pappert.  It  was  used  at  a  distance 
of  4  km  from  the  transmitter. 
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IV.  EXAMPLES 


Three  examples  are  given  of  the  implicit-finite-difference  method.  The  two  models  that  were  success¬ 
fully  used  were  each  of  a  surface  duct. 

The  first  model  was  of  a  bilinear  profile,  such  that  the  upper  part  had  a  slope  of  dM/dz  =  0.118  and 
the  lower  part  had  a  slope  of  dM/dz  =  -0.118.  The  apex  was  at  14  m,  the  transmitter  was  at  25  m,  and 
the  receiver  at  5  m.  Starting  values  were  found  at  4  km  by  using  the  ray-trace  formulation.  The  vertical 
mesh  size  was  0.25  m  and  the  horizontal  step  size  was  100  m.  An  absorber  was  used  from  96  through  128 
m.  Path  loss  as  a  function  of  distance  is  shown  in  figure  1  as  the  solid  line.  The  x  symbols  in  the  figure 
are  waveguide  normal  mode  theory  values  as  found  in  the  program  MLAYER.  It  is  seen  that  the  values 
agree  very  closely. 

The  second  example  is  of  a  multilinear  profile  representing  a  14-m  surface  duct  [from  Hitney,  1988]. 
The  profile  is  given  in  table  1.  The  transmitter  was  at  25  m  and  the  receiver  at  5  m.  The  “free-space” 
starter  was  used  at  2.5  km  from  the  transmitter.  The  vertical  mesh  size  was  0.25  m  and  the  horizontal  step 
size  was  20  m.  An  absorber  was  used  from  75  through  100  m.  Path  loss  as  a  function  of  distance  is  shown 
in  figure  2.  The  x  symbols  are  waveguide  normal  mode  values.  It  is  seen  that  the  values  agree  very 
closely. 


0.00  9.00  18.00  27.00  36.00  45.00  54.00 

RANGE  (km) 

Figure  1.  Bilinear  model. 
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V.  SUMMARY  AND  CONCLUSIONS 

The  theory  of  the  implicit-finite-difference  method  has  been  summarized  and  the  boundary  conditions 
described.  Examples  indicate  that,  for  sufficiently  high  signal  levels,  the  IFD  method  is  adequate  at  9.6 
G Hz.  However,  the  method  did  not  give  correct  results  beyond  35  km  for  the  standard  atmosphere  case  in 
which  signal  levels  are  very  low. 
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